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Abstract. The "loop equations" of random matrix theory are a hierarchy of equations born of attempts to 
obtain explicit formulae for generating functions of map enumeration problems. These equations, originating 
in the physics of 2-dimensional quantum gravity, have lacked mathematical justification. The goal of this 
paper is to provide a complete and short proof, relying on a recently established complete asymptotic 
expansion for the random matrix theory partition function. 



1. Background and Preliminaries 

The study of the Unitary Ensembles (UE) of random matrices begins with a family of probability 
measures on the space of N x N Hermitian matrices. The measures are of the form 

dut = J- cxp {-N Tr [V t (M)]}dM, 
Zn 

where the function Vt is a scalar function, referred to as the potential of the external field, or simply the 
"external field" for short. Typically it is taken to be a polynomial, and written as follows: 



1 - 

V t = V(X; t 1 ,...,t v ) = -\ 2 + ^t j \ j 



2 

where the parameters {t\, . . . , t v } are assumed to be such that the integral converges. For example, one may 
suppose that v is even, and t v > 0. The partition function Zn, is the normalization factor which makes the 
UE measures be probability measures. 

Expectations of conjugation invariant matrix random variables with respect to these measures can be 
reduced, via the Weyl integration formula, to an integration against a symmetric density over the eigenvalues 
which has a form proportional to below: 



(1.1) exp^-TV 



1 N 1 
-]TV(A 3 ; h,...,t v )- —J2log\\j-\ e 



d N X. 



These latter multiple integrals can be more compactly expressed in terms of kernels constructed from polyno- 
mials {pi (A)} orthogonal with respect to the exponential weight e~ NVt ^ For instance, the fundamental 
matrix moments E (TrAPJ, where E denotes expectation with respect the measure d/it, are expressed as 

/oo 
\ j (\)d\ 
-oo 

where p^ (A) denotes the so-called one-point function 
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p%\x) = (i# {j : A,- e (-00,1)} 



(1.3) 



1 

iV 



K N (x,x) 



with the kernel Kjsr{x,y) 



e 2 



(Vt(s)+Vt(2/)) 



^p £ (x)^(y). 



The symbol E denotes expectation with respect to the normalization of the measure (|1.1|) which is given by 
dividing this family of measures by the corresponding family of integrals: 



(1.4) 



Z]\r(ti,t2, ■ ■ ■ , t v ) 



exp < —N 



1 N 1 



d N \. 



We will sometimes refer to the following set of t = (ti, . . . , t v ) for which l|1.4|l converges. For any given T > 
and 7O, define 

v-l 

T(r,7)-{tel": |t|<T, * u >7X)l*il}. 

The leading order behavior of Z/v(ii, t<z, . . . ,t v ) is rather classical, and is known for a very wide class of 
external fields V (see, for example, ^O]). We will require the following result. 

Theorem 1.1. There is T > and 7 > so that for all t £ T(T, 7), the following holds true: 

(1) 



(1.5) 

where 
(1.6) I(t u ...,t v ) 



J™ ■^2 1 °g{^V-(*li*2 J • ■ -,tv)} = ~I(h, ...,t v ) 
N— >oo iv 



inf 



V(A)d/x(A) - / / log|A-/,|4KA)d/*fo) 



Borel measures p,p>0,/ dfj,— l 

(2) There is a unique measure fj,y which achieves the infimum defined on the right hand side of 111.6$ . 
This measure is absolutely continuous with respect to Lebesgue measure, and 

dfxv = ip d\, 

V-(A) = ^X(^)(A) v / (A-a)(/3-A)MA), 

where h(X) is a polynomial of degree v — 2, which is strictly positive on the interval [a, (3\ (recall that 
the external field V is a polynomial of degree v ). The polynomial h is defined by 

h(z) - 1 i *"« dS 



2tti J v /(i - a)y/{s - 0) s- z 

where the integral is taken on a circle containing (a, [3) and z in the interior, oriented counter- 
clockwise. 

(3) There exists a constant I, depending on V such that the following variational equations are satisfied 
by Hv ■ 



(1.7) 



J 2 log |A - n^dnvW) + V{\) > I for A e R\supp(fi v ) 
2 log |A — 77| -1 G^iy(?7) + V(A) = I for Ae supp(/j, v )- 



(4) The endpoints a and [3 are determined by the equations 



f 



V'(s) 



ds = 







V(s-a)(/3-s) 
sV'(s) 



ds = 2tt. 



y/(8-a)(J3-8) 



(5) The endpoints a(t) and /3(t) are actually analytic functions oft, which possess smooth extensions 
to the closure of {t : t e T(T, 7)}. TTiey afeo satisfy — a(0) = /3(0) = 2. /n addition, the coefficients 
of the polynomial h(X) are also analytic functions of t, with smooth extensions to the closure of 



Remark The variational problem appearing in ({1.6(1 is a fundamental component in the theory of random 
matrices, as well as integrable systems and approximation theory. It is well known, (see, for example, . 
that under general assumptions on V, the infimum is achieved at a unique measure /iy, called the equilibrium 
measure. For external fields V that are analytic in a neighborhood of the real axis, and with sufficient growth 
at 00, the equilibrium measure is supported on finitely many intervals, with density that is analytic on the 
interior of each interval, behaving at worst like a square root at each endpoint, (see 0] and □ 

Remark For a proof of l|1.5[) , we refer the reader to ^U] > however this result is commonly known in the 
approximation theory literature. □ 

Remark It will prove useful to adapt the following alternative presentation for the function if>: 



where the function R(X) is defined viai?(A) 2 = (A — a) (A — f3), with R(X) analytic in C\ [a, /?], and normalized 
so that -R(A) ~ A as A — > 00. The subscript ± in i?±(A) denotes the boundary value obtained from the upper 
(lower) half plane. □ 

The goal of this paper is to provide a rigorous justification for the "loop equations" which originated in 
the physics of 2-dimensional quantum gravity (see, for example, the survey |d] and the references contained 
therein). More precisely, this entails 

• Proving that the quantity Pw _ dz possesses a complete asymptotic expansion in even powers 
of N (see Theorem IO and 

• Establishing that the coefficients P g (z) in the asymptotic expansion l|3.16fl satisfy the hierarchy of 
nonlocal equations (|3.18l) . which are the loop equations. 

Once the coefficients P g {z) are known to exist as analytic functions of z and the times t, they may be 
interpreted as generating functions for a collection of graphical enumeration problems for labelled maps, 
counted according to vertex valences and the genus of the underlying Riemann surface into which the maps 
are embedded. (See and also [?]•) Because of the combinatorial connection and its use in 2-dimensional 
quantum gravity, obtaining explicit formulae has been a fundamental goal within the physics community of 
quantum gravity. The loop equations arose as a means to obtain explicit information (and possibly explicit 
formulae) for these coefficients, although without mathematical justification. 

In Section |3 we will need to consider the Cauchy transform of the equilibrium measure: 



{t : t 6 T(T,7)}, with 



h(X,t = 0) = 1. 



(1.8) 



1&(A) = — R+(X)h(X), AG (a, /J), 

ZTTl 




z e C/R. 
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It follows from differentiating the variational equations, Theorem ll.lf 3'). that F(z) solves the scalar Riemann- 
Hilbert problem 

F + (s)+F_(s) = V'(s), s€[a,/3] 
F + (s)-F_(s) = 0, seR/[a,0\ 

F(z) = - z +0{z- 2 ), z^oo. 

From this it is straightforward to deduce that 

(1.9) 2F(z) = V'(z)-^-R(z)h(z) 

2. Large N Expansions 

The fundamental theorem for establishing complete large N expansions of expectations of random variables 
related to eigenvalue statistics was developed in [J]. A concise statement of this result is: 

Theorem 2.1. [Ercolani and McLaughlin, [7J/ There is T > and 7 > so that for all t G T(T, 7), the 
following expansion holds true: 

/oo 
f(X)p$(X)dX = fo + N- 2 h + J\r 4 / 2 + • • • , 
-00 

provided the function /(A) is C°° smooth, and grows no faster than a polynomial for A — » 00. The coefficients 
fj depend analytically on t fort G T(T, 7), and the asymptotic expansion may be differentiated term by term. 

The complete details for the derivation of this result are presented in [7J; however, there are a few specifics 
presented there that we repeat here for use in subsequent sections and for general background information: 

• The function p^ has a full and uniform asymptotic expansion, which starts off as follows: 

(2.2) /9 W(A)=^(A) + o(7V- 1 /2). 

• The specific form that this expansion takes depends very much on where one is looking; for example, 
for A G (a, (3), the expansion takes: 

(2.3) = ip(X) + JL (_L_ _ _L_) cos Lf m4 



1 

+ N2 





H(X) + G(X) sin ^ N I ip(s)ds 



A 



in which H (A) and G(A) are locally analytic functions which are explicitly computable in terms of 
the original external field V(X). 
In [7] the primary application of this theorem was to establish a complete large N symptotic expansion 
of II. II exists: 

Theorem 2.2. There is T > and 7 > so that for t G T(T, 7), one has the N — > 00 asymptotic expansion 
(2.4) log = N 2 e (x,t) + ei (x,t) + ^e 2 (x,t) + ■■■. 

The meaning of this expansion is: if you keep terms up to order N~ 2h , the error term is bounded by CN~ 2h ~ 2 , 
where the constant C is independent t for all t G T(T, 7). For each j, the function ej(t) is an analytic 
function of the (complex) vector (t), in a neighborhood of (0). Moreover, the asymptotic expansion of 
derivatives oflog(Zjv) may be calculated via term-by-term differentiation of the above series. □ 

Remark Recently, Bleher and Its [3] have carried out a similar asymptotic expansion of the partition 
function for a 1-parameter family of external fields. A very interesting aspect of their work is that they 
establish the nature of the asymptotic expansion of the partition through a critical phase transition. □ 



Remark A subsequent application in [Sj is to develop a hierarchy of ordinary differential equations whose 
solutions determine recursively the coefficients e g for potentials of the form V(X) — A 2 /2 + tX 2v . 



Remark The asymptotic results in |7| were also used recently in [S] to establish that asympotics of each 
individual eigenvalue have Guassian fluctuations, regardless of whether one is in the bulk or near the edge 
of the spectrum (provided only that the eigenvalue number, when counted from the edge, grows to oo). 



In the present paper we will make use of a mild extension of Theorem 12. II in which the function / is of 

w{ 
X- 



the form with z living outside of the inteval [a, [3\ 



Theorem 2.3. For each 5 > 0, there is T > and 7 > so that for all t € T(T, 7), the following expansion 
holds true: 

(2.5) y°° (y^j P^(A)dA = w (z) + N^w^z) + N- 4 w 2 (z) + ■■■ , 

provided z € C \ [a — 6, [3 + S], the function w(X) is analytic in a neighborhood of R, and grows no faster than 
a polynomial for X — > 00. The coefficients Wj depend analytically on z and t for t S T(T, 7), and possess 
convergent Laurent expansions for z — > 00. Furthermore, the asymptotic expansion may be differentiated 
term by term. 

For z bounded away from the real axis, this Theorem follows from Theorem 12 .11 The mild extension to 
the case when z may be near the axis (but bounded away from the support [a, (3\) follows by exploiting 
analyticity to replace the integral along the real axis near z by a semi-circular contour so that A remains 
uniformly bounded away from z. Once the contour is such that A is bounded away from z, the uniform 
asymptotic expansion for may be used. Since z is away from the support [a, 0\, the newly introduced 
semi-circular contour is also bounded away from the support, and one may use arguments similiar to those 
presented in \J\ Observation 4.2] (where they were used for A real and outside the interval [a,/3]) to show 
that (A) is uniformly exponentially small on the semi-circular contour, and also that the residue term 
obtained from deforming the contour is also uniformly exponentially small. We will leave these details for 
the interested reader. 

3. Derivation of the Loop Equations 
We introduce some notation. Denote the Greens function of a random matrix M as 

G(z,M) = (z-M)- 1 

and and its trace as 

g{z) = Tr G(z). 
We evaluate dG^i / dMij in two different ways to get a useful relation: 
Lemma 3.1. 

E(GkiGji) — N E(GkiV (M)ji) 

Proof. Since G ■ (z — M) — 1, d/dM^ G ■ (z - M) = 0; or, equivalently 

dG/dMij ■ (z - M) - G ■ E i3 = 0, 

where Eij is the elementary permutation matrix with a 1 in the ij th entry and all other entries zero. It 
follows that 

dG/dMij = G ■ E i:j ■ G; 

in particular, 

dG H /dMij = {G ki Gji), 

and so 

E (dGki/dMij) = E (GkiGji) . 
5 



On the other hand, integrating by parts yields 

EidGu/dMij) = J- / ^exp{-NTr[V t (M)]}dM 
Z N Jn oMij 

= -J- / G kl7 ^-exp{-NTr[V t (M)]}dM 
Z N Jn oM i:j 

= I G k iTr(VV(M)-E ij )exp{-NTr[V t {M)]}dM 

Z N Jn 

= N ±- f G k iV'(M) 3l exp{~NTr[V t (M)}}dM 
Z N Jn 

= NE (GkiV'(M)ji) . 
Combining the above two representations for E (dGu/dMij) gives the result. □ 
Proposition 3.2. 

(3.1) E((g(z)f)=NE(Tr(G-V'(M))) 

This follows directly from the lemma by setting i — k and j = I, summing over k and I and dividing by 
1/N 2 which yields 

J2^(GkkG u ) =Nj2^(Gk,iV'(M) lk ) 

k,l k,l 



The relation (|3.1() can be naturally regarded as a generating function for the second order matrix cumulants 
of M when written in the equivalent form 

(3.2) E((. 9 (z)) 2 ) - E((. 9 (z))) 2 = AE (±Tr (G ■ V'(M)^j E((. 9 (z))) 2 . 

To proceed further we will need to introduce some more notation. First, we will use the following general 
expression for the potential 

oo 
3=0 

which is understood to have only finitely many tj non-zero. We also have the formal vertex operator 
(3 ' 3) dV = ~^^+idL- 

3=0 3 

(A precise meaning for this formal relation will be given in the beginning of the next section.) This can be 
used to give a compact formal representation of a generating function for matrix moments in terms of the 
RM partition function 

^ wh** z » = £^ E (^ TrM 0- 

3=0 V 7 

3.1. Asymptotic Expansions. In order to make formal relations such as (|3.4|l meaningful we need to use 
some fundamental asymptotic facts. The trace of G(z) has a standard integral representation in terms of 
the RM one-point function p. 3(1 

(3.5) g(z) = N T ^-d\. 



z-X 

By boundedness and exponential decay of (A), g(z) has a valid asymptotic expansion in large z as 



(3.6, ~ e> >«o 

j=0 



Thus (|3.4[) can be precisely understood as saying that for each m and large z 

i , „ ^ i r°° av^(a) 



d y(m) N 2 & ^+1 / Z - A 



dX 



^^dA + 0(z-< m+1 )) 
z — A 



where 

d 



m—1 ^ 



d y(m) z j+i dt 

J=0 1 



In a similar sense we have the following asymptotic equation for each m and large z: 

m 1 

E((.g(z)) 2 ) - E((.g(z))) 2 = ]T (e (w • TrM fc ) - E (TrAP) E (TrM fc )) + 0(«" 



j,fe=0 



1 -logZ N + 0(z-^+^) 



dV(m) dV( m ) N 2 

(3. 8) . ;j JL T £^ +0(2 -(- + 3» ) . 

In what follows, we will use instead of d y\ m) but with the above asymptotic interpretation understood. In 
the rest of this section we need to establish that there are estimates controlling the errors in the asymptotic 
expansions (|3.7(l and l|3.8|) that remain valid uniformly as N — > oo. To this end we first note that for l|3.7|) 
the error has the form 



i [°° A">^(A) dA 



z m+1 z-X 



±# {/<">(*) + N^f[ m \z) + N-\ft\z) + ■■■} 



for t £ T(T, 7). The RHS is a uniformly (in T) valid asymptotic expansion which follows from the funda- 
mental Theorem 12. II . Similarly for l|3.8(l the error has the form 

00 1 A 1 

- E ^T^^T {ftH*) + N-*f[ m \z) + N-±ft\ Z ) + ■■■} 

k—m 

E {^M + ^ 2 <?i m) W + at 4 ^*) + • • • } 



in which the sum on the RHS is finite since V depends on only finitely many distinct . We use here the fact 
stated in Theorem 12 . II that these asymptotic expansions in TV can be differentiated term by term, preserving 
uniformity. 

With the observations of this section we may express the relation (|3.2|l in terms of integral representations 
involving the one-point function: 



dVJ.^ z-X J_ OQ z-X \J-00 z-X J 

to be understood in the sense of an asymptotic expansion in large z whose coefficients moreover have 
asymptotic expansions in even powers of iV which are uniform in admissible t. We note that one consequence 
of this is that the two terms on the RHS of (|3.9|l cancel at leading order so that the difference has leading 
order 0(1) for large N. 

3.2. Loop Equations. To prepare for the transformation to a recursive loop equation we parse the first 
integral on the RHS of (|3.9II as 

(3 . 10) f!Mffl fl . rnMM d , +0{ e-°»), 

J -00 z-X J a _ s z-X 

7 



where c > depends on the choice of the positive constant 6. The justification for the exponential error 
term is part of the proof of the fundamental asymptotic relation presented in f. ■ 
From (|3.1U|) .we may further transform this term: 



z - X 2ni J a ^ s \J C x - X J z-X 

(3.12) 



2ni J c J a _s \ x — X J z — \ 

0+5 



1 

(3.13) = ^££V VV g) W (_l_)(_l_) (iA<tt + ( e -* ) 

2tti Jc z- x J a _ s x — X 

where in the first line we have expressed V'(X) as a loop integral a la Cauchy's Theorem, with the contour 
C encircles the interval (a — 6, (3 + 6), with z outside the contour of integration, and in 1)3. 13)1 . one term has 
vanished by Cauchy's Theorem and analyticity. 

Inserting 1)3. 14|) into 13.9fl . we have derived the final form of the Loop Equation generating function. 



z — x J a —$ x — A dV j a _ ^ z — A \ j a _5 z — A 



Using Theorem 12 .31 the term J^f P " Ja ^ dX is easily seen to possess an asymptotic expansion in large N, 
each of whose coefficients possesses a Laurent expansion in large z: 

(3.16) / esMax^N-^p^). 

Ja-6 g=Q 

Combining (|3.7() and (|2.4|) we see that 

. . d , . d , , 1 de n (t) 

(3-17) PM - ^) = ^r 1 )^(t) = -E i iTr^ i - 

j=0 J 

Inserting 1)3.16)1 into the Loop Equation generating function 1)3.15)1 yields the hierarchy of Loop Equations: 

(3.18) ±.J^ Pg{x)dx = _^_p 5 _ l(z )_^p 5 ,( z) p g _ 3 , (z) 

m Jc z x g , =0 

JC g'=i 

where the transition to the final recursion formula is mediated by the identity ()1.9)) 

(3.19) ip(x) = V'(x) - 2 r ^-dX = V'(x) - 2P (x) 

J -co X — X 

where ip( x ) here is interpreted as the analytic extension of the density for the equilibrium measure off of the 
slit [a,0\ as given by 1)1.8)1 . 



With this result in hand it is now possible to consider a recursive derivation of the terms P g starting 
with Pn as given by (|3.19() . These terms are related to the map enumeration generating functions, e s (t), 
through l|3.17fl . I n the physics literature there are instances in which loop equations are used to formally 
derive expressions for some of the e g . In particular, we refer the reader to pQ. 

A natural application of our derivation of (|3.18[) would be to the derivation of closed form expressions for 
e g (t) which extend our results in [5] for potentials V depending only on a single non-zero time ti v . We may 
also be able to use these equations to say something about the general qualitative and asymptotic behavior 
of the e g . Finally, the P g (z,t) contain information about the large N asymptotic behavior of the matrix 
moments, such as (|1.2|) . which could be used to explore whether or not the general unitary ensembles are 
asymptotically free. 
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